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Abstract 

Fermat's principle and variational analysis is used to analyze the trajectories of light prop- 
agating in a radially inhomogeneous medium with a singularity in the center. It is found that 
the light trajectories are similar to those around a black hole, in the sense that there exists a 
critical radius within which the light cannot escape, but spirals into the singularity. 

1 Introduction 

There has recently, within the scientific community, been much interest focused on the extra- 
ordinary properties associated with Bose-Einstein condensates - clouds of atoms cooled down to 
nano-Kelvin temperatures where all atoms are in the same quantum state and macroscopic quan- 
tum conditions prevail ||l|, ^, |^ . 

One of the most dramatic experiments which has been performed in Bose-Einstein condensates 
is the demonstration of optical light pulses traveling at extremely small group velocities, e.g. 
velocities as small as 17 m/s have been reported |^. 

In a physically suggestive application of the optical properties of Bose-Einstein condensates, it 
was recently suggested that they could be used to create, in the laboratory, dielectric analogues of 
relativistic astronomical phenomena like those associated with a black hole |Q. 

The astrophysical concept of a black hole is one of the most fascinating and intriguing phe- 
nomena related to the interaction between light and matter through the curvature of space-time 
within the framework of general relativity. Black holes are believed to form when compact stars 
undergo complete gravitational collapse due to, e.g., accretion of matter from its surroundings (in 
general relativity, the pressure contributes to the gravitational mass of a fluid, and an increased 
pressure will, after a certain point, therefore only help to accelerate the collapse phase). 

In ordinary dielectric media, unrealistic physical conditions would be required in order to 
demonstrate any of the spectacular effects of general relativity. However, in dielectric media 
characterized by very small light velocities, new possibilities appear. In particular, it was recently 
suggested that by creating a vortex structure in such a dielectric medium it is possible to mimic 
the properties of an optical black hole |H . What this study in fact suggested was that one could 
construct, by using the aforementioned vortex, an unstable photon orbit, much like the orbit at 
a radius r — 3M in the Schwarschild geometry |^. In order to construct an event horizon, one 
would need to supplement the vortex flow by a radial motion of the fluid |||, |^ . 

The vortex, which involves a rotating cylindrical velocity fleld, tends to attract light propagating 
perpendicular to its axis of rotation and to make it deviate from its straight path. In the model 
considered in Refs. |^, the increasing, in fact diverging, rotational velocity field of the vortex 
core attracts light by the optical Aharanov-Bohm effect and causes a bending of the light ray 
similar to that of the gravitational fleld in general relativity. It is also shown that there exist a 
critical radius. Verity from the vortex core with the properties that if light rays come closer to the 



core than Tcrit, the Hght will fall towards the singularity of the vortex core. In this sense, the critical 
radius, rcrit, plays the role of an optical unstable photon orbit, analogous to the unstable orbit in 
the Schwarzschild geometry. 

The basic physical effect involved in the analysis of the ray propagation in Refs. ||l|, |^ is the 
fact that the refractive index, n, of a medium changes when the medium is moving. In Refs. |0, ^, 
the medium is assumed to have a cyHndrical vortex velocity field V — V(r) given by 

W 

V(r) = -0 (1) 

where r is the radius from the vortex core centre, ip is the azimuthal unit vector and 2ttW is the 
vorticity. 

Advanced physical concepts like Bose-Einstein condensates, vortex velocity fields and black 
holes are fascinating physical effects, which however may be difficult to present in a simplified 
manner for undergraduate students. Nevertheless, it is important to try to convey inspiration for 
and arouse curiosity in such phenomena. An interesting example of such an eflFort was recently 
made in Ref. where the basic physical mechanism, in the form of a classical two level system 
was used to demonstrate the possibility of the very low light velocities observed in Bose- Einstein 
condensates. In a later work by the same author, Ref. |^, an inspiring investigation is given of a 
simple mechanical model that exhibits a gravitational critical radius. The purpose of the present 
work is similar to that of Refs. ^ ^. We will try to discuss a simple classical example of a medium 
where the refractive index has a divergence at r = and which exhibits some of the dramatic 
properties of the light behaviour around a black hole as discussed in Refs. ||l|, H, ||, ^. In addition 
we want to illustrate the power and beauty of the classical principle of Fermat and of variational 
methods in connection with this new fascinating concepts. 



2 Fermat 's principle 



It is well known that Hght tends to be defiected towards regions with higher refractive index. 
Consider light propagating in a cylindrically symmetric medium where the refractive index increases 
towards the centre. In such a situation we expect the light path to look qualitatively as in Fig. 1. 
The actual light path is determined by Fermat's principle i.e. 

(2) 



S J n{r)ds = 

whereds is an infinitesimal element along the light ray. We will consider light propagation in a 
medium where the refractive index is of the qualitative form, cf Eq. (1) 



n(r) 



1 as r ^ ro 
ro/r as r <C ro 



One possible realization of n{r) with the desired properties is given by 



l+(^ 
r 



Using Fermat's principle and expressing the light path as the relation 9 
polar angle, Eq. (2) implies 

S J n(r)^l + r2 (^^^ dr = . 
The Euler-Lagrange variational equation corresponding to Eq. (5) reads 

de/dr 



d_ 

dr 



n{r)- 



+ [de/dr) 
which determines the trajectory of the light. 







(3) 
(4) 

(r) where is the 
(5) 

(6) 
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3 Solution of the light trajectory 

Equation (6) directly implies that 

r'^n(r) dO/dr 
^1 + r2(rf0/rfr)2 



(7) 



where the constant is determined by initial conditions, i.e. the properties of the incident hght. 
Equation (7) is easily inverted to read 

— = ± (81 
dr ^r^n^{r)-rfr^ 

For a light ray incident as shown in Fig. 1, we clearly have d6/dr < (at least up to some 
minimum radius r = r„ii„. It is illustrative to first consider the trivial case of a homogeneous 
medium with n{r) = 1. In this case Eq. (9) becomes 

dr - rfr^ ^ ' 

which can easily be interpreted to yield the hght path in the form 



6* = arccoty ^ - 1 (10) 

or simpler 

r- = ^ (11) 

Clearly this is the straight line solution 

y = ri (12) 

where the parameter plays the role of "impact parameter" or minimum distance from the centre. 

Let us now consider the model variation for n(r) as given by Eq. (5). Within this model Eq. 
(10) becomes 

- = (13) 
dr ^r^ - (r? - r§)r2 

Clearly the solution of Eq. (13) will depend crucialh^ on the relative magnitude of r-i and tq, i.e. 
the impact parameter relative to the characteristic radial extension of the inhomogeneity. Consider 
first the case when tq < r^. Equation (13) can then be rewritten as 



6i de 6i 



n dr i/r4 - 62^2 



(14) 



b, ^ (15) 
Equation (15) is of the same form as Eq. (9) and we directly infer the following solution 



sm 



(16) 



As ^ — » 0, we still have asymptotically 



r ~ ^ (17) 



However, the trajectory is now bending towards the origin and the minimum distance occurs at 
the polar angle 6 = 6^ given by 
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The corresponding minimum distance, is 



- . /^2 



] (19) 

We also note that the trajectory is symmetrical around the angle 0^ and that the asymptotic 
angle of the outgoing light ray is 



lim 9{r) 



26„ 



(20) 



The solution given by eq.(lO) describes a trajectory which is bent towards the centre of at- 
traction at r = 0. Depending on the ratio r^jri, the trajectory is either more or less bent or may 
even perform a number of spirals towards the centre before again turning outwards and escaping, 
cf. Fig. 2. The number of turns, A^, which the light ray does around the origin before escaping is 
simply 



N 



10 m 




1 


_ 2tt _ 







(21) 



where [xj denotes the largest integer less than x. 

Let us now consider the special case when the impact parameter equals the characteristic width 
of the refractive index core i.e. = tq. The equation for the trajectory now simplifies to 



d9 

dr 



with the simple solution 



n 

r 



(22) 



(23) 



i.e. the trajectory describes a path in the form of Arkimede's spiral as the light falls towards the 
origin, cf. Fig. 3. The form of the light trajectory in the situation when < ro is now obvious, 
it will spiral into the singularity more or less directly depending on the magnitude of the ratio 
ro/n > 1. 

The actual trajectory in this case is given by a slight generalization of Eq. (16) viz 



sinh[0-\/l - rf/r^] 



(24) 



This solution does indeed convey the expected behaviour, the trajectory spirals monotonously 
into the singularity of the refractive index. 



4 Final comments 

The present analysis is inspired by recent discoveries and discussions about light propagation in 
Bose-Einstein condensates where extremely low light velocities can be obtained. This has triggered 
speculations about possible laboratory demonstrations of effects, which normally are associated 
with general relativistic conditions. In particular, it has been suggested that it may be possible to 
create the analogue of a black hole using a divergent in-spiral of a Bose-Einstein condensate. 

In the present work we have analyzed a simple classical example of light propagation as deter- 
mined by Fermat's principle in a medium characterized by a radially symmetric refractive index. 
In analogy with the variation of the vortex velocity field suggested in |^ , the refractive index is here 
assumed to diverge towards the centre. This classical example exhibit some of the characteristic 
properties of light propagating around a black hole where the gravitational attraction deflects the 
light and where, under certain conditions, the light may be "swallowed" by the black hole. 

In the example analyzed here, the unstable photon orbit of the Schwarzschild black hole is 
similar to the characteristic radius, tq, of the refractive index variation, which together with the 
impact parameter, r^, of the incident light completely determines the light trajectory. If tq < Vi, 
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the light is more or less deflected, but ultimately escapes. However, if ro > Vi, the light spirals 
into the singularity. It should be cautioned that this result depends crucially on the presence of 
the singularity in the refractive index (4). If this is removed the in-spiraling photon orbit will 
eventually turn and start spiraling outwards. In the recent discussion about the possibility of 
generating analogues of optical black holes in Bose-Einstein condensates, it has been suggested 
that the vortex motion must also have a velocity component in the radial direction. 

In the classical model considered here, a number of additional physical efi^ects will obviously 
affect the light path when it comes close to the axis and will in fact remove the mathematical 
singularity. Nevertheless, the model provides a simple example of light dynamics, which resembles 
some of the properties of light propagation around a black hole. Furthermore, the investigation is 
based on Fermat's principle and variational analysis, in this way illustrating the use of classical 
methods in connection with very new and fascinating concepts at the front Hne of modern research. 

Figure Captions. 

Fig.l Qualitative plot of a light ray trajectory in a cylindrically symmetric medium with a refrac- 
tive index, which increases towards the centre. 

Fig. 2 Light trajectories for the refractive index model of eq.() for different impact radii > ro 
Fig. 3 Light trajectory in the case of Tj = ro, the spiral of Arkimede. 
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